Measurement of Stochastic Entropy Production 
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Using fluorescence spectroscopy we directly measure entropy production of a single two-level 
system realized experimentally as an optically driven defect center in diamond. We exploit a recent 
suggestion to define entropy on the level of a single stochastic trajectory (Seifert, Phys. Rev. Lett. 
95, 040602 (2005)). Entropy production can then be split into one of the system itself and one of 
the surrounding medium. We demonstrate that the total entropy production obeys various exact 
relations for finite time trajectories. 



Entropy as the central concept in statistical physics 
pervades many branches of science. Well-defined and 
uncontested only in equilibrium, its extension to time- 
dependent non-equilibrium phenomena has been debated 
since the days of Boltzmann mostly in relation to an ex- 
planation of irreversibilityand a foundation of the second 
law of thermodynamics |]|. Major progress arose with 
the formulation of the fluctuation theorem, which quan- 
tifies in the long time limit the probability of entropy an- 
nihilating trajectories in small systems constantly driven 
in a steady state 0, 0, 0, 0, @ • Entropy production in 
these systems is either defined as phase space contraction 
rate or associated with a dissipation functional, which ul- 
timately should describe the dissipated heat. 

By introducing the notion of a stochastic entropy along 
a single trajectory, it has become possible both to ex- 
tend the validity of the fluctuation theorem to finite times 
and to prove an integral fluctuation theorem for the to- 
tal entropy production in arbitrarily driven systems gov- 
erned by stochastic dynamics 0]. In this Letter, using 
our previously introduced driven two- level system 0] , we 
measure the stochastic entropy production along single 
trajectories and demonstrate that it obeys various exact 
relations for finite times. 

Fluctuation theorems for entropy production should be 
distinguished from related theorems like the Jarzynski re- 
lation 
lation 



Crooks' theorem [lfj , and the Hatano-Sasa re- 
11| . The first two allow to extract free energy dif- 
ferences from non-equilibrium work measurements. Ex- 
perimental tests have been performed by using a tor- 
sional pendulum ^3], by mechanically stretching RNA 
hairpins [lH H3 | as well as by driving a colloidal particle 
in a time-dependent harmonic [l5lllq | and non-harmonic 
potential |l7]]. The third one yields an exact relation for 
transition between different steady states from which a 
general Clausius inequality follows that has been tested 
using a driven colloidal particle All these relations 
address a small system embedded in a surrounding heat 
bath of constant temperature. In contrast, our set-up 
works athermally and does therefore neither involve nor 
require any notion of dissipated heat. While our previ- 
ous work has demonstrated an exact Jarzynski-like rela- 
tion for the athermal analog of "dissipated work" 0] , the 



present paper addresses the concept of stochastic entropy 
production directly. The crucial difference is that the 
derivation of the latter requires using the actual nonequi- 
librium probability distribution in a general master for- 
mula 0|, whereas the former involves the corresponding 
equilibrium distribution fl9j . Still, both quantities fulfill 
various exact relations like an integral fluctuation theo- 
rem. 

Our system is a photochromic defect center in natural 
Ila-type diamond. Its optical properties indicate that we 
are dealing with a nickel-related center j2jj. It can be 
excited by red light responding with a Stokes-shifted flu- 
orescence. Additionally to this "bright" state the defect 
exhibits a non-fluorescent "dark" state. The transition 
rates a (from dark to bright) and b (from bright to dark) 
depend linearly on the intensity of green and red light, 
respectively, turning the defect center into an effective 
two-level system 



(dark) ^ 1 (bright) 

b 



(i) 



with controllable transition rates a and b. 

Single centers were addressed with a home-built con- 
focal microscope [2l| using a dye laser (CR699, DCM) 
running at 680 nm (red light intensity oc b) superim- 
posed to the 514 nm line of an Ar-Ion laser (green light 
intensity oc a) as excitation sources. While the red light 
was kept constant throughout the experiment, the green 
light was modulated using a function generator controlled 
acoustooptical modulator ( AOM) . In addition to the flu- 
orescence of the single defect, a second Avalanche Photo 
Diode recorded simultaneously the alternating intensity 
of the green light. 

The system can be found in state n with probability 
p n , where n takes cither the value or 1. To drive the 
system out of equilibrium, we modulate the rate a (from 
dark to bright) according to the sinusoidal protocol 



a(t) = a [l +7sin(27rt/t m )] 



(2) 



whereas the rate b is held constant. The parameters are 
the equilibrium rates ao and b, the period t m , and the 
modulation depth ^ 7 < 1 . At a time resolution of 1 ms 
the data of the two detectors were acquired simultane- 
ously after starting the modulation protocol for the green 
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FIG. 1: Entropy production in a single two-level system with parameters ao = (15.6ms) _1 , b — (21.8 ms) _ , t m = 50ms, and 
7 = 0.46. (a) Transition rate a(t) (solid/green line) and probability of the bright state pi(t) (red circles - measured, dashed/red 
line - theoretical prediction) over 4 periods, (b) Single trajectory n(t). (c) Evolution of the system entropy. The gray lines 
correspond to jumps (vertical dotted lines) of the system whereas the dark lines show the continuous evolution due to the 
driving, (d) Entropy change of the medium, where only jumps contribute. (e,f) Examples of (e) entropy producing and (f) 
entropy annihilating trajectories. The change of system entropy As = s(t) — s(0) (solid/black) fluctuates around zero without 
effective entropy production, whereas in (e) As m (dashed/red) produces a net entropy over time. In (f), As m consumes an 
entropy of about 1 after 20 periods, (g-i) Histograms taken from 2,000 trajectories of the system (g), medium (h), and total 
entropy change (i). The system entropy shows four peaks corresponding to four possibilities for the trajectory to start and end 
(0 h 1, 1 h 0, h 0, and 1 i— » 1). The distribution (h) of the medium entropy change has mean (As m ) = 1.7 and width 
<7 m = 3.7, on this scale it differs only slightly from the distribution of the total entropy change (i). 



laser. Following 20 periods with t m = 50 ms and a certain 
modulation depth 7, the system was given 1,000 ms of 
unmodulated green light to relax back into equilibrium. 
A sine function has been fitted to the mean intensity of 
the green laser to obtain the modulation depth 7. The 4 
runs in this Letter contain 3 times 1,000 and once 2,000 
trajectories. The modulation depths were varied from 
7 = 0.07 to 0.46. 

A dimensionless, non-equilibrium entropy for driven 
systems on the level of a single stochastic trajectory has 
been defined in Ref. 01 as 

»(t) = -ln[p n(t) (f)] , (3) 

where the measured probability p n is evaluated at the 
actual state n(t) at time t. Fig. shows the protocol 
a(t) together with the probability pi(t) to dwell in the 
bright state. Fig. ^> displays a sample binary trajec- 
tory n(t) jumping between the two states. In Fig. ^ we 
see that the evolution of s(t) is governed by two effects. 
First, the time-dependent driving of the rates leads to an 
evolving probability resulting in a continuous contribu- 
tion. Second, a jump between the two states gives rise 
to a contribution — ln[p+/p_], where p_ and p+ are the 
probabilities of the states immediately before and after 
the jump, respectively. 



Beside the entropy of the system itself, energy ex- 
change and dissipation lead, in general, to a change in 
medium entropy. For an athermal system, this change in 
medium entropy As m can not be inferred from the ex- 
changed heat. Rather it has to be defined. In Ref. 0, 
the choice 

A Sm = ln^ (4) 

Wji 

for a jump from state i to state j with instantaneous 
rate tuy (Wji being the backward rate) has been mo- 
tivated in analogy to the thermal case. In our case it 
becomes As m = — ln[a(i)/6] for a jump 1 1— ► and 
As m = —ln[b/a(t)] for a jump 1— ► 1. As demon- 
strated in Fig. ^i, the medium entropy changes only 
when the system jumps, thereby balancing to some de- 
gree the change of s(t). 

One of the fundamental consequences of a stochastic 
entropy (0 is the fact that beside entropy producing tra- 
jectories also entropy annihilating trajectories exist, see 
Fig. and f, respectively. However, in accordance with 
physical intuition, the latter become less likely for longer 
trajectories or increased system size. Entropy annihilat- 
ing trajectories not only exist, they are essential to satisfy 
the integral fluctuation theorem [j| 

(exp[-As t ot]) = 1. (5) 



3 





FIG. 3: Experimental test of the detailed fluctuation theo- 
rem (DFT) ® after 20 periods with the same parameters as 
in Fig0 Following Eq. (|SJ, the slope should be 1. The in- 
set shows the slope of the DFT as function of the trajectory 
length. Due to the statistics of the entropy distribution, the 
DFT is satisfied best for longer trajectories. 
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FIG. 2: Experimental test of the integral fluctuation theo- 
rem JSJ. (a) Probability distribution of the total entropy pro- 
duction of a driven two-level system. The modulation depth 
7 increases from 0.07, 0.14, 0.20 to 0.46; other parameters as 
in Fig.0 The trajectory length is 20 periods with 1,000 tra- 
jectories per distribution (2,000 trajectories for 7 = 0.46). (b) 
The mean (exp[— Astot]) over trajectory length for the four 
different modulation depths. 



This theorem states that the average (• • •) over infinitely 
many realizations of a process involving the total entropy 
change As to t = As + As m becomes unity for any trajec- 
tory length and any driving. Trajectories with As to t < 
may occur seldom but are exponentially weighted and 
thus contribute substantially to the left hand side of 
Eq. JSJ. As an immediate consequence of Eq. J5J one 
has with (As to t) ^ a consistent formulation of the sec- 
ond law of thermodynamic for small systems, giving a 
posteriori support to the entropy definition ■ 

Fig. EK demonstrates how the average entropy produc- 
tion increases with increasing driving amplitude 7. To 
fulfill the constraint imposed by Eq. © the distributions 
spread, making trajectories with large negative produc- 
tion more likely. In Fig. |2Jd we present the experimen- 
tal evidence for the validity of the theorem JSJ analyzing 
1,000 trajectories. Only for the largest modulation depth 
in combination with long trajectories (t > 15 periods) a 
deviation from the theorem is observable. This is due 
to the need for la rger statistics as the mean value of the 
entropy increases [23j |. 

A stronger but also more special version of a fluctu- 
ation theorem is the detailed fluctuation theorem 



P(+As tQt ) 
P(-Aa to t) 



exp[+As to t]- 



(6) 



Its experimental test is shown in Fig. [3J Adapted to our 
situation, it states that entropy annihilating trajectories 
with probability P{— As to t) become exponentially less 
probable compared to trajectories producing the same 
amount of entropy with the increasing absolute value of 
the entropy. This theorem was derived originally for the 
long-time limit in non-equilibrium steady states. How- 
ever, following Crooks' reasoning 01 it even holds for 
periodic driving as in our experiment, provided (i) the 
protocol is time-symmetric, and (ii) the distribution p n {t) 
has relaxed into the corresponding periodic limit distri- 
bution. The first requirement is easily met by choosing 
start and end of the trajectories at n + 1/4 periods. As 
for the second requirement, in Fig. ^1 we see that the 
probability p n (t) indeed relaxes and trails with a con- 
stant phase shift behind the driving rate a(t). Therefore 
we wait several periods before we start to record the data. 

In conclusion, we have provided the first direct mea- 
surement of entropy production along single stochastic 
trajectories in periodically driven systems. In particular, 
we have shown that it is crucial to include the entropy of 
the system itself. While for long trajectories it remains 
bounded, its contribution is required for both the inte- 
gral and the detailed fluctuation theorem to hold for fi- 
nite times. How this stochastic entropy will contribute to 
a more comprehensive understanding of non-equilibrium 
dynamics remains to be elucidated through the future 
study of more complex systems. 
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